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ABSTRACT 


In  recent  years,  the  development  of  mathematical  models 
for  the  study  of  ground  shock  effects  in  soil  and/or  rock 
media  has  made  important  progress.  Currently,  three  basic 
types  of  advanced  models  have  been  studied:  (a)  elastic- 

ideally  plastic  models,  (b)  variable  moduli  models,  and 
(c)  elastic- non ideally  plastic  capped  models.  This  paper 
is  the  first  part  of  an  investigation  which  is  aimed  at 
studying  the  behavior  of  the  different  types  of  models  in 
ground  shock  problems. 

The  ground  shock  response  in  the  superseismic  range  of 
a  one  megaton  air  burst  on  a  homogeneous  half-space  of  a 
soil  is  considered.  Each  of  the  three  tvpes  of  models  was 
fitted  to  laboratory  test  data  and  calculations  were  made 
for  each  case.  Generally,  the  results  from  all  three  models 
are  comparable  only  when  the  stress  paths  in  uniaxial  strain 
are  comparable  for  complete  load-unload  cycles.  When  this  is 
not  so,  major  differences  occur  in  the  lateral  motions  and 
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INTRODUCTION . 


In  recent  years,  the  development  of  mathematical 
models  for  the  study  of  ground  shock  effects  in  soil  and/or 
rock  media  has  made  important  progress,  Refs.  [1],  [2]. 

Currently,  three  basic  types  of  advanced  models  have  been 
studied:  (a)  elas t  ic-ideally  plastic  models;  (b)  variable 

moduli  models;  and  (c)  elastic-nonideally  plastic  capped 
models.  A  detailed  summary  of  the  advantages  and  disadvantages 
of  each  type  of  model  is  given  in  Refs.  [1],  [21.  Table  I  of 
Refs,  [1],  [2]  is  repeated  in  this  report  for  convenience. 

The  present  paper  is  the  first  part  of  an  investigation 
aimed  at  studying  the  behavior  of  the  different  types  of 
models  in  ground  shock  problems  of  interest.  *he  study  in¬ 
volves  fitting  a  mathematical  model  of  each  type  to  the  same 
set  of  experimental  data  and  then  using  these  models  to 
analyze  a  series  of  ground  shock  problems.  The  aim  of  the 
investigation  is  to  evaluate  the  influence  of  the  various 
constitutive  models  on  ground  motion  calculations,  and  to 
explain  any  differences  that  might  arise.  It  is  intended  to 
expand  this  study  to  include  comparative  calculations  of 
field  tests  and  eventually  to  make  recommendations  on  the  use 
of  the  models  for  a  vide  range  of  pressure  levels,  geometries 
(layering)  and  soil  properties. 

This  paper  considers  a  one  megaton  (1  Ml')  air  burst  on  a 
homogeneous  half-space,  Fig.  (1),  of  a  soil  material,  namely 


McCormick  Ranch  sand.  This  material  was  ctosen  because  it  was 
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one  of  the  very  few  soil  materials  for  which  a  large  amount 
of  laboratory  data  was  available.  Material  model  fits  for 
the  two  more  sophisticated  results,  i.e.,  the  variable  moduli 
model,  Ref.  [3],  and  the  soil  cap  model,  Ref.  [4],  were  already 
available  from  our  previous  studies.  The  elas t ic-ideally 
plastic  fit  was  developed  specifically  for  the  present  study. 

Results  are  monitored  for  a  number  of  surface  pressure 
levels  in  the  high  superseismic  region  (Mach  =  8)  and  a 
comparison  of  the  ground  shock  results  obtained  with  all  three 
models  is  made.  As  might  be  expected,  the  results  indicate 
the  great  importance  of  modeling  correctly  the  stress  path  in 
unloading  in  uniaxial  strain.  It  should  be  noted  that  in  the 
past,  stress  path  information  from  the  uniaxial  strain  test 
was  usually  missing  from  available  experimental  results.  One 
purpose  of  this  paper  will  be  to  emphasize  the  importance  of 
making  such  information  available  to  the  model  developer  by 
means  of  a  uniaxial  strain  null  test.  Important  studies 
towards  obtaining  reliable  null  test  data  are  currently  in 
progress  at  the  Waterways  Experiment  Station. 

The  soil  selected  for  modeling  was  McCormick  Ranch  sand 
because  this  was  one  of  the  few  materials  for  which  almost 
all  of  the  laboratory  data  from  the  standard  tests  was 
available.  This  data  included:  (a)  uniaxial  strain  test 

results  for  Longitudinal  stress  versus  longitudinal  strain 
over  a  load-unload  cycle  (however,  no  lateral  stress  versus 
longitudinal  stress  data  and  hence,  no  stress  path  was 
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available) ;  (b)  triaxial  compression  test  results  including 

stress  difference  versus  strain  and  strain  difference  curves  1 

.  1  •  j  ,  ( 

in  cyclic  loading  and  stress  difference  at1  failure  for  each 

test;  and  (c)  proportional  loading  test  axial  stress  versus 

1  1  ^  f 

axial  strain  results.'  Th$  materia!  is  characterized  by  a 

considerable  amount  of  hysteresis  in  a  ‘  1 oad^unload  cycle  in 

either  dniaxial  strain  or  triaxial. compression . 

■  * 

The  study  presented  here  is  limited  t'**  the  superseismic 

region  of  the  materia].  Generally-  the  results  from  all  three 

■ 

models  are  comparable  when  the  stress  pa<ths'in  uniaxial  strain 

{  1 

are  comparable,  for  complete  load-unlo!ad  cycles.  ^hen  this  is 
not  so,  major  differences  can  be  found  in  the  lateral  motions  f 
and-  stresses,  although1  the,  vertical  mot'ions  and  stresses  are 
still  governed  primarily  by  the  stress-strain  curve  for  the 
uniaxial  strain  test  and  hence,  are  only  slightly  affected. 
Consequently ,.  particular  attention  i.,  paid  ts  the  evaluation 
of  these  lateral  (horizontal)  motions  in  this  pap>er. 

>  , 

This  first  study  does  not  ‘include  layering  effe'ets.  Some 
preliminary  results  for  the  case  of  atrigid  bottom  at  a  given' 

i 

depth  are  also  discussed.  The  study  is  currently  being 

i 

extended  to  other  geometries  (layered  materials) ,  range^  » 

(  transseismic  and  subseismic)  and  materials  of  interest. 


i 
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II  MATERIAL  MODELS. 


McCormick  Ranch  sand  was  used  previously  to  illustrate 

i 

the  fitting  of  both  the  variable  moduli  and  soil  cap  models 

i 

to  a  real  material.  Each  of  these  fits  was  done  by  different 
investigators,  independently ,  using  the  same  rather  extensive 

i 

set  of  laboratory  data.  The  descriptions  of  these  fits, 
including  comparisons  with  experimental  data,  are  well 
documented  in  previous  reports,  Ref.  [3],  for  the  variable 
moduli,  and  Ref.  [4]  for  the  soil  cap  model.  The  minor 
differences  in  the  fits  to  ,  t  h  e  various  test  curves  are  not 
significant.  The  important  thing  to  note  is  that  each  model 
represented  a  good  fit  to  the  available  laboratory  data. 

Many  of  the  differences  that  occur  are  attributable  to 

i 

differences  in  procedure  and  philosophy  of  the  individual 
modeler.  For  example,  the  constant  unloading  (elastic)  bulk 
modulus  in.  the  cap  model  is  not  necessary  to  the  cap  model 
formulation.  A  variable  unloading  bulk  modulus,  such  as  is 
used  in  the  variable  moduli  approach,  could  easily  be 

f 

incorporated  into  the  cap  model.  The  elas tic-idea lly  plastic 
fit  was  made  specifically  for  this  comparison.  It  cannot 
really  be  considered  to  be  completely  independent  of  the 
variable  moduli  fit,  since  it  was  made  by  the  same  person 
after  the  variable  moduli  fit  had  already  been  completed. 

In  fact,  certain  relations  such  as  the  unloading  bulk  modulus 
•  * ) 

and  the  failure  envelope  are  identical  in  these  two  models. 


The  "failure  envelope"  for  the  variable  moduli  material 
can  be  defined  as  the  surface  f(J^  ,  p)  =  0  such  that 

gld  =  °- 


f 
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McCormick  Ranch  sand  is  characterized  by  a  large  amount 
of  hysteresis  in  a  load-unload  cycle  in  uniaxial  strain,  see 
Fig.  (2).  Also  illustrated  in  Fig.  (2)  are  stress-strain 
relations  in  uniaxial  strain  computed  with  the  three  models. 
Exce;  t  for  the  tail  at  low  stress  in  unloading  both  the 
variable  moduli  and  cap  fits  follow  the  experimental  curve 
extremely  well.  The  elas t ic-ideally  plastic  fit  which  was 
used  deviates  somewhat  more.  The  reader  is  referred  to 
Refs.  [3]  and  [4]  (or  more  briefly  to  Sections  IV  and  V  of 
Refs.  [1],  [2])  for  additional  laboratory  data  as  well  as 

the  details  of  the  fits.  A  brief  description  of  the  three 
types  of  models  follows. 

( A )  Elas t ic-Ideally  Plastic  Model . 

Model  of  this  type  are  currently  being  used  for  most 
ground  shock  calculations.  They  have  been  described  in 
detail  in  several  papers,  Refs.  [1],  [2],  [5],  [6],  and  only 

a  summary  is  given  here. 


The  model  is  characterized  by  a  yield  condition  of  the 
type,  Fig.  (3) 


k  +  3ap(l 


=  k  +  f  apc 


eff 


P  <  Pc 
P  >  Pc 


(1) 


For  this  soil,  a  nonassociat ed  flow  rule  of  a  von  Mises  type 
is  used  to  remove  "dilatancy"  effects  from  yielding.  This 
leads  to  a  unique  pressure-volumetric  strain  relation  on 


initial  loading.  This  pressure,  which  is  given  by 


^max^kk^  ~  an  £kk 
n  =  l 


is  the  maximum  pressure  possible  at  a  given  volumetric 
strain.  For  unloading,  and  subsequent  reloading  up  to  the 
previously  reached  maximum  pressure  on  the  virgin  loading 
curve,  a  variable  bulk  modulus  of  the  form 


(2) 


Kuidp^  "  Kou  +  Kiup 


=  K 


max 


P  iP. 


P  -  Pm 


(3) 


is  used.  In  the  model,  Poisson's  ratio  is  assumed  to  be 
piecewise  constant,  i.e., 


V  =  =  Constant  (Loading) 

(4) 

V  -  =  Constant  (Unloading  and  Reloading) 

where  the  choice  is  based  upon  which  pressure-volume 
relation  is  used.  In  the  calculations,  a  tension  cutoff 
was  used 


p  >_  0  (No  mean  tension)  (5) 

Two  models  were  used  in  the  calculations:  (1)  one  in  which 
the  same  constant  value  of  V  was  used  in  both  loading  and 
unloading  and  (2)  one  in  which  different  constant  values 
of  and  VyN  were  used  for  loading  and  unloading,  respective¬ 

ly.  A  summary  of  the  material  constants  used  in  the  elastic- 
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*) 

ideally  plastic  fit  is  given  in  Table  II. 

As  is  usual  in  such  models,  the  stress-strain  curve  in 
uniaxial  strain  was  modeled  over  a  full  load-unload  cycle 
and  the  yield  condition  was  fitted  to  the  failure  envelope 
from  the  triaxial  compression  tests.  The  stress  difference- 
strain  curves  for  the  triaxial  tests  could  not  be  modeled 
once  yielding  started;  this  is  a  limitation  on  and  a 
characteristic  of  such  elas ti c-idea lly  plastic  models. 

(B )  Variable  Moduli  Model. 

Models  of  this  type  have  been  discussed  in  detail  in 
Refs.  [1],  [2],  [3],  [7],  [8],  and  som^  two  dimensional 
calculations  have  been  made,  Ref.  (9],  for  comparison  with 
a  series  of  analytical  solutions,  Refs.  [10],  [11].  In  this 
model,  both  the  bulk  and  shear  moduli  are  functions  of  the 
stress  and/or  strain  invariants  and  no  explicit  yield 
condition  is  specified.  The  model  gives  good  qualitative 
agreement  with  the  usually  available  soil  laboratory  test 
data.  The  mathematical  description  is  given  in  terms  of  the 
incremental  stress-strain  relations 


p  =  K  e 


kk 


s .  .  =  2G  e ,  , 
ij 


(a) 

(b) 


(6) 


It  should  be  noted  that  the  separation  of  the  constitutive 
relations  into  a  deviatoric  and  a  volumetric  part,  while 


*) 


It  should  be  noted  that  not  all  of  the  constants  in 
Table  II  are  independent. 
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convenient,  precludes  dilatancy  in  the  material.  Conse¬ 
quently,  this  model  is  not  intended  to  represent  materials 
such  as  rock  when  significant  dilatancy  may  be  observed. 

The  pressure-volume  relations  are  given  by 

3 

pmax (ekk}  =  l.  \  Ekk  Initial  loading 

n  =  l 


kun(p)  K0U  +  Klyp 


=  K 


p  >  0 


max 


Tension  Cutoff 


p  <  p 
r  —  rm 


P  -  Pm 


(7) 


(8) 


(9) 


The  quantity  K^Cp)  is  used  for  both  unloading  and 
reloading  regimes. 


The  deviatoric  stress-strain  behavior  is  given  by 
Eq.  (6b)  where  the  shear  modulus  G  is  defined  separately 
in  loading,  unloading  and  reloading  as 


GLD(P>VJ2) 


G  =  GUN(P*  VJ2) 


GR£(p>  VJ2> 


when  F  =  F  and  J0  >  0 
max  2 


when  £  0 


when  F  <  F  and  J0  >  0 
max  2 


(a)  N 

(b) 

(c)  ) 


>  (10) 


where  the  quantities  G^  ,  G^  are  given  by 


LD  K  ’ 

=  G1 
I  t 

JUN  ^P ’  V J2^  =  G0U  1U1 

=  G1U 


2Ul 


P  <  P, 


P  >  P, 


P  1  Pc 


>  (ID 


p  >  p 


c  / 


During  reloading  (in  shear)  defined  by  Eq,  (10c),  the  shear 


! 

f 


modulus  G__  is  taken  as  the  linear  combination  of  GT n  and 


RE 


LD 


UN 


GRE^P,~VJ2^  ^F  ^  GLD<P>  VV  *  ^  ~  ^F  ^  ^  GUN^P,~VJ?^ 

max  max 


2 ' 

(12) 


where  n  is  a  constant.  The  ratio 


max 


is  a  measure  of  the 


current  state  of  stress  F,  defined  by 


t  '-'t  n 

F(P,  -\/J,)  =  1  -  - 


G,n<P»  V^) 


GLD(P’  0) 


(13) 


to  the  maximum  previous  state,  F  .  The  function  F  is 
r  max 


restricted  to  the  range 


0  <  F  <  1  (14) 

so  that 

gld(p*  Vv  -  GRE(p>  -  gun(p>  (15) 

The  various  constants  }  used  in  this  fit  are  given  in 

*  *  ) 

Table  III.  Detailed  comparisons  of  this  fit  to  laboratory 
data  for  McCormick  Ranch  sand  are  presented  in  Ref,  [3]  and 
in  Section  IV  of  Refs.  [1],  [2]. 


J  As  in  the  previous  case,  not  all  of  these  constants  are 
independent . 

£  jSf  \ 

J  The  current  variable  moduli  fit  is  a  slight  modification 
of  the  fit  referred  to  as  "Uniax-Triax  II"  in  Ref.  [3]. 
The  unloading  shear  parameters  Gq^  ,  (and  y£jj)  are 

somewhat  altered.  The  reloading  parameter  n,  not  given 
elsewhere,  was  chosen  as  4  somewhat  arbitrarily. 
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The  theory  of  variable  moduli  models  is  rigorously 
correct  for  cases  of  proportional  loading  in  shear.  Con¬ 
tinuity  problems  may  arise  for  complicated  stress  paths  in 
which  neutral  loading  in  shear  occurs.  One  purpose  of  this 
study  is  to  ascertain  the  behavior  of  the  model  in  ground 
shock  problems  which  although  not  a  case  of  proportional 
loading  in  shear,  would  not  be  expected  to  exhibit  significant 
neutral  loading. 

(C)  Soil  Cap  Model. 

This  model  represents  the  latest  step  in  the  development 
of  soil  and  rock  models,  Refs.  [1],  [2],  [4],  [12],  It  is 

characterized  by  the  specification  of  elastic -non ideally 
plastic  constitutive  relations  with  a  capped  yield  surface. 

A  yield  surface  which  combines  both  ideal  plasticity  and 
strain-hardening  is  postulated  and  an  associated  plastic  flow 
rule  is  used.  The  model  is  capable  of  predicting  observed 
laboratory  data  from  all  available  tests  and  at  the  same  time 
satisfies  the  most  restrictive  theoretical  requirements  of 
continuity,  stability  and  uniqueness. 

For  soil,  the  yield  conditions,  Fig.  (4),  consist  of  an 

*  ) 

ultimate  failure  condition  of  a  form  slightly  different 
than  Eqs.  (1) 

fl(Jl  ’  V5  =  VJ2  "  [A  "  C  exP(-BJl>]  =  0  (16) 

where  J  =  3p,  and  an  interval  strain-hardening  cap  which 

'  This  difference  is  another  example  of  difference  in 
philosophy  of  the  two  modelers  referred  to  earlier. 

Equations  (1)  could  have  been  used  here  just  as  well. 
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"expands"  or  "contracts"  as  the  plastic  volumetric  strain 
increases  or  decreases,  respectively.  In  the  present  fit  to 
McCormick  Ranch  sand,  an  elliptical  cap,  See  Fig*  (4),  of  the 
form 

f2(Jl  »  eP)  =  (J1  "  L)2  +  ~  (x  -  L)2  =  0  (17) 

was  used,  where  R  is  the  constant  ratio  of  the  major  to  the 
minor  axis,  and  where  L(e  )  and  X(e  )  are,  respectively,  the 
values  of  at  the  center  of  the  ellipse  and  at  the  inter- 

p 

section  of  the  cap  with  the  axis.  The  function  Xt,e  ) 
relating  the  position  of  the  cap  to  the  plastic  strain  is 
given  by 

X(eP)  •=  -  [£n  (1  -  eP/W)]/D  (18) 

where  W  and  D  are  constants* 

Within  the  current  yield  surface  the  material  behaves 
elastically  with  a  constant  shear  modulus  G.  The  elastic 
bulk  modulus  K  can  be  a  function  of  pressure  without 
raising  any  theoretical  objections.  However,  for  the 
present  fit  to  McCormick  Ranch  sand  a  constant  bulk  modulus 
was  considered  sufficient.  Consequently,  the  value  of 
Poisson's  ratio  within  the  yield  surface  also  is  a  constant. 

The  values  of  the  constants  used  in  the  current  fit  are 
given  in  Table  IV. 

All  theoretical  requirements  for  continuity,  stability 
and  uniqueness  are  satisfied  by  the  cap  model.  It  should 
Ke  noted  that  since  the  movement  of  the  cap  is  controlled 
by  the  increase  or  decrease  in  the  plastic  volumetric  strain, 
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No  attempts  were  made  to  make  the  models  similar  and 
hence,  each  fit  represents  x^hat  a  qualified  investigator 
might  come  up  with  independently,  given  the  laboratory  data 


available.  This 

aspect  of 

the 

model  fitting 

will  be 

quite 

important  in  the 

discussion 

of 

the  numerical 

results 

in 

Section  III. 


The  use  of  a  pressure  dependent  yield  condition  with  an 
associated  flow  rule  necessarily  leads  to  "dilatancy " , 
an  inelastic  volume  change  due  to  yielding. 
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III  COMPUTATIONAL  RESULTS  AND  DISCUSSION. 

For  comparative  purposes,  the  models  were  each  separately 
run,  using  the  LAYER  Code,  for  the  case  of  a  1  MT ,  1500  ft 
height  of  burst  explosion  on  a  hair-space  of  the  McCormick 
Ranch  sand  material.  In  the  computations,  the  surface  loading 
(Brode  H-l  loading,  Ref.  [13]),  the  number  of  grid  points  in 
the  horizontal  and  vertical  directions,  the  horizontal  and 
vertical  spacing  of  the  grid  points  (20  ft  by  20  ft)  and  the 
treatment  of  the  grid  boundaries  were  identical  for  all  three 
runs.  The  left  boundary  included  the  origin,  R  =  0,  the  right 
boundary  assumed  the  velocities  to  be  continuous  across  the 
boundary  while  the  bottom  boundary  was  of  the  transmitting 
type,  Refs.  [14],  [15],  so  that  a  half-space  wat  effectively 
simulated.  All  runs  used  an  iterative  stress  boundary 
condition  on  the  surface,  z  =  0,  Ref.  [15],  so  that  the 
applied  surface  stress  and  the  constitutive  relations  were 
both  satisfied.  Of  course,  different  constitutive  relations 
were  applied  for  each  run. 

Since  the  unloading  bulk  modulus  in  both  the  variable 

moduli  and  the  elas tic-ideally  plastic  models  varied  with 
pressure,  the  maximum  unloading  P  wave  velocities  at  high 
pressures  (4157  ft/sec  and  4368  ft/sec,  respectively)  were 
considerably  greater  than  the  constant  (elastic)  unloading 
P  wave  speed  of  2068  ft/sec  for  the  cap  model.  Consequently, 
while  the  run  for  the  cap  model  could  be  made  with  6  ms  time 
steps,  Courant-Friedr ichs-Lewy  stability  required  the  time 
step  for  the  other  two  runs  to  be  only  3  ms. 
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For  the  highly  sup er seismic  range  being  studied  here, 
the  vertical  displacements  were  considerably  larger  than  the 
horizontal  displacements.  Detailed  results  are  presented  for 
the  displacements  over  a  substantial  portion  of  the  super- 
seismic  region  and  in  particular  at  the  1000  psi  (R  =  1660  ft) 
and  500  psi  (R  -  2060  ft)  surface  pressure  levels. 

Figure  (5)  shows  a  plot  of  the  horizontal  displacement 
versus  range  curves  for  all  three  models  at  time  t  =  1167  ms 
for  a  near  surface  point  (depth  z  =  30  ft).  It  is  seen  that 
while  the  horizontal  displacements  from  the  variable  moduli 
and  cap  models  are  comparable  (within  20%),  the  elastic- 
ideally  plastic  model  gives  significantly  ^mailer  horizontal 
displacements,  sometimes  by  a  factor  of  as  much  as  two  (when 
compared  with  the  variable  moduli  results).  On  the  other  hand, 
a  similar  plot  of  vertical  displacement  versus  range  gives 
comparative  results  which  are  virtually  indistinguishable 
for  all  three  models. 

To  investigate  the  reason  for  the  differences  in  the 
highly  sensitive  horizontal  displacements,  the  time  histories 
of  the  horizontal  velocities  at  the  500  psi  contour  are  shown 
in  Fig.  (6)  for  the  cap  and  the  elas tic-ideally  plastic  models. 
It  is  seen  that  while  at  early  times  these  histories  are 
quite  comparable,  the  signal  from  the  cap  model  computations 
remains  at  a  considerably  higher  level  than  that  from  the 
elas ti c-ideally  plastic  model  computations  at  later  times; 
this  results  in  a  significantly  larger  horizontal  displacemen t 
for  the  cap  model.  To  further  understand  these  differences, 
the  time  plots  for  the  radial  stress  0^  are  shown  in  Fig.  (7) 
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for  both  models.  Again  it  is  seen  that  for  very*  early  times 

1  i 

the  signals  are  comparable  but  that  at  later  times  significantly 
larger  radial  stresses  are  given  by  the  cap  model. 

The  reason  for  these  differences  is  best  seen  by  considering 
the  stress  paths  of  the  three  models  in  a  uniaxial  strain  con¬ 
figuration  for  a  complete  load-unload  cycle.  For  the  present 

highly  superseismic  case  (Mach  number  -  8)  the  actual  state 

<  ■  » 

of  stress  in  the  material  can  be  considered  similar  to  that  ini 

,  i 

a  uniaxial  strain  test. 

,  '  / 
i 

The  stress  path  in  uniaxial  strain,  where  there  is  a 

*  *  t 

single  independent  stress  deviator,  is  illustrated  in  Fig.  (8) 

as  a  plot  of  the  axial  stress  deviator  s  versus  the  mean  stress 

or  pressure  p.  In  this  (s^  s  p)  space,  points  on  a  line  with  i 

a  slope  of  -1  correspond  to  a  constant  value  of  0^  ,  while 

those  on  a  line  with  a  slope , of  +2  correspond  to  a  constant 

value  of  0^  .  The  loci  0^  «  0  and  0^  «  0  are  shown  in  Fig.  (8). 

Points  on  lines  parallel  to  a  =  0  (a  =  0)  but  to  the  right 

z  r 

correspond  to  la,rger  values  of  a  (a  )  .  The  failure  envelope, 

z  r 

f  , 

£qs.  (l)or  Eq,  (16),  appears  on  the  plot  as  a  pair  of  curves 

rr-  /j 

symmetric  with  respept  to  the  p , axis  since  y 3^  -  +  ~  for 
this  geometry. 

f  f  : 

The  stress  paths  for  all  three  models  s>tart  at  the  origin 
(stress  free  state)  and  load  up  to  an  axial  stress  0 =  400  psi. 

The  cap  and  the  variable  moduli  curves  are  very  close  to  one 

.  { 

another,  starting  almost  tAngent  to  0  =0  (which  wpuld 

j  r  ' 

correspond  to  Vq  =  0) ,  and  bending  over  as  they  approach  the 
failure  surface.  The  elas tic-ideally  plastic  model,  with  a 


r 
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constant  Poisson’s  ratio,  loads  along  the  straight  line  \2J 
which  falls  short  of  intersecting  the  failure  envelope.  The 
unloading  path  for  this  model  differs  even  more  drastically 
from  those  of  the  other  two  models. 

The  curve  marked  (?)  is  also  the  unloading  stress  path 
for  the  original  1  elas  tic-ideally  plastic  model  x^hich  utilized 
the  same  constant  Poisson’s  ratio  V  =  0.30  in  both  loading 
and  unloading.  It  is  seen  that  the  material  loads  along  a 
straight  line  to  G^  =  400  psi  and  unloads  back  down  the  same 
line  to  G^  =  0.  At  the  completion  of  unloading,  the  radial 
stress  G^  is  equal  to  zero. 

The  dotted  curve  shows  the  stress  path  in  loading  and 
unloading  for  the  cap  model.  The'  loading  path  is  a  curve  up 
to  the  400  psi  level,  and  a  straight  line  for  unloading  to 

f  i 

i 

the^  zero  level  at  the  point  B.  It  is  important  to  note, 

however,  that  the  cap  model  (which  is  linearly  elastic  in 

ualoading)  unloads  with  a  constant  value  of  Poisson’s  ratio 

\)tt„  =  0.25  and  that  when  the  a  =0  line  is  reached  at 
UN  z 

point  B,  a  residual  radial  stress  a  of  75  psi  is  maintained 
in  the  material. 

The  variable  ihoduli  model,  represented  by  the  thick 

f 

black  line  in  the  figure,  loads  along  a  curve  to  G  =  400  psi 
and  then  unloads  along  a  different  curve  to  the  u  =0  line 
at  the  point  C  where  a  level  of  residual  radial  stress  of 
G  =44  psi  is  noted.  It  is  now  apparent  as  to  why  the  radial 
stress  from  the  ground  shock  calculations  for  the  cap  and  the 
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elas tic-ideally  plastic  models  differed,  since  the  former 
unloads  quite  differently  from  the  latter  (i.e.,  it  maintains 
a  residual  radial  stress).  This  difference  in  the  stress 
path  for  the  three  models,  and  in  particular,  the  major 
differences  in  the  unloading  stress  path  Setween  the  cap 
and  the  e las t ic-j dea lly  plastic  models  are  responsible  for 
the  large  differences  in  the  horizontal  displacements  which 
were  shown  in  Fig.  (5)  at  the  500  psi  pressure  contour.  The 
loading  stress  paths  and  more  significantly,  the  unloading 
stress  paths  for  the  cap  and  the  variable  moduli  model  are 
quite  similar,  thus  leading  to  radial  stresses  and  velocities 
which  are  similar  and  hence,  more  comparable  horizontal 
displacements,  as  shown  in  Fig.  (5).  The  elas tic-ideally 
plastic  model  has  a  quite  different  stress  path  history  and 
hence,  gives  significantly  different  (in  this  case  lower) 
radial  stresses,  velocities  and  displacements. 

To  test  this  marked  dependence  of  the  horizontal  quantities 
on  the  stress  paths  in  loading,  and  more  particularly  unloading, 
a  second  calculation  was  run  with  an  elas t ic-i deally  plastic 
model  in  which  the  Poisson's  ratio  in  unloading  was  matched 
with  that  of  the  cap  model,  i.e.,  =  0.25.  The  loading 

stress  path  for  this  case  is  again  given  by  the  straight  line  (?) 
in  Fig.  (8),  but  the  unloading  curve  is  now  parallel  to  that  of 
the  cap  model  and  is  shown  by  the  straight  line  (T) .  It  is 
noted  that  upon  complete  unloading  to  =  0,  a  residual  radial 
stress  =  38  psi  [Point  A  in  Fig.  (8)]  is  obtained.  It  will 
be  shown  later  in  this  section  that  the  horizontal  displacements 
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for  this  revised  elas tic-ideally  plastic  model  (in  whi-h  the 
unloading  stress  path  is  more  comparable  to  those  for  the  cap 
and  the  variable  moduli  models)  are  much  closer  to  the  corre¬ 
sponding  values  for  the  cap  and  the  variable  moduli  models. 

The  importance  of  having  models  which  give  similar  stress 
paths  in  loading-unloading  cycles  is  further  illustrated  by 
considering  a  loading  to  the  higher  stress  level  O »  1000  psi, 
as  shown  in  Fig.  (9).  In  this  case,  the  loading  was  sufficiently 
high  for  the  elas tic-ideally  plastic  model  to  plasticize.  The 
elas tic-ideally  plastic  model  loads  along  the  straight  line, 
yields  at  the  intersection  of  this  line  with  the  yield  surface 
and  continues  loading  along  the  yield  surface  to  the  a  ^  =  1000  psi 
line.  It  then  unloads  linearly  along  line  2  (when  =  0.30 

is  used)  or  line  1  (when  =  0.25  is  used)  until  it 

replasticizes  on  the  other  side  of  the  yield  surface  and  then 
continues  along  the  lower  yield  surface  until  the  final  point 
at  Oz  =  0  is  reached.  The  cap  model  loads  along  the  dotted 
curve  to  the  =  1000  psi  line,  and  then  unloads  along  the 
straight  dotted  line  until  the  line  meets  the  lower  portion  of 
the  yield  surface  and  replasticization  occurs.  It  then  unloads 
along  the  lower  yield  surface  to  the  final  point  at  =  0. 

The  variable  moduli  model  loads  along  the  heavy  black  curve 
and  unloads  along  the  heavy  black  curve  until  once  again  the 
final  point  at  =  0  is  reached. 

For  this  case  of  loading  to  a  high  stress  level,  i.e.,  a 
level  in  which  the  yield  surface  plays  a  significant  role  in 
the  elas ti c-ideally  plastic  model,  all  three  models  have  a  very 


19 


similar  loading  stress  path.  However,  the  stress  paths  in 
unloading  are  still  different,  but  not  nearly  as  different  as 
they  were  for  the  previous  case,  Fig.  (8).  One  might  therefore 
expect  that  the  computed  horizontal  displacements  at  the  1000  psi 
level  for  the  elas tic-ideally  plastic  model  with  =  0.3,  the 

cap  model  and  the  variable  moduli  model  would  agree  somewhat 
better  than  those  for  the  case  of  loading  to  500  psi.  This  is 
in  fact  borne  out  in  Fig.  (5).  For  the  case  of  the  elastic- 
ideally  plastic  model  with  =  0.25  (equal  to  for  the 

cap  model),  the  stress  path  histories  are  quite  similar  for 
all  three  models  and  much  closer  agreement  would  be  expected 
for  the  horizontal  displacements  from  the  ground  shock 
calculations,  as  shown  in  the  curves  which  follow. 

Figures  (10)-(20)  show  comparative  results  for  the 

calculations  with  the  elas tic-ideally  plastic  model  with 

VtTXT  =  0.25,  the  cap  model  and  the  variable  moduli  model.  It 
UN 

is  seen  that  generally,  excellent  agreement  is  obtained  from 
the  three  models.  Figure  (10)  shows  the  curve  of  vertical 
displacement  versus  range  for  all  three  models.  It  is  seen 
that  the  results  are  practically  indistinguishable.  If  the 
uniaxial  strain  test  stress-strain  relations  for  the  three 
models  were  identical  over  the  entire  stress  range  shown  in 
Fig.  (10),  the  three  curves  in  the  figure  would  essentiaJ ly 
coincide.  While  the  curves  produced  by  the  variable  moduli 
and  elas tic-ideally  plastic  models  almost  coincide,  the 
vertical  displacements  computed  with  the  cap  model  are  slightly 
larger  at  both  the  high  and  low  stress  levels.  The  three  curves 
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do  coincide  at  about  the  1000  psi  overpressure  level  since  the 

models  were  all  fit  using  the  experimental  uniaxial  strain 

test  to  1000  psi,  Fig.  (2).  No  attempt  was  made  to  make  the 

models  do  more  than  behave  reasonably  when  extrapolated  to 

stress  levels  significantly  above  1000  psi,  since  uniaxial 

strain  data  was  not  available.  The  unloading  bulk  moduli  for 

the  variable  moduli  and  e las tic-ideally  plastic  models  are 
*) 

identical  ,  while  for  peak  vertical  stresses  less  than  about 
370  psi  (p  =  240  psi)  the  cap  unloading  bulk  modulus  is  stiffer 
than  that  of  the  other  models  and  consequently,  the  vertical 
displacement  is  larger. 

Figure  (11)  shows  the  variation  with  range  of  the  maximum 
horizontal  displacements  computed  with  the  three  models.  It 
is  seen  that  the  peak  horizontal  displacements  are  fairly  close 
over  the  entire  range,  with  the  maximum  difference  being  on 
the  order  of  20%.  The  results  for  the  elas tic-ideally  plastic 
model  are  consistently  lower  than  those  for  the  other  two 
models.  The  variable  moduli  and  cap  results  agree  extremely 
well  over  a  wide  range  of  surface  overpressures. 

Figure  (12)  shows  the  attenuation  of  peak  vertical  and 
radial  stresses  with  depth  below  the  (nominally)  300  psi 
surface  pressure  contour.  The  computed  vertical  stress 
attenuations  for  the  three  models  agree  extremely  well  for 
depths  less  than  z  -  150  feet.  Below  that  level,  the  additional 
hysteresis  in  the  cap  model,  due  to  the  stiffer  unloading  bulk 
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These  two  models  were  fit  by  the  same  investigator. 


modulus  at  low  stresses,  causes  the  additional  attenuation 
in  the  vertical  stress  for  that  model  to  be  noticeable.  Near 
the  surface,  where  the  peak  vertical  stresses  agree,  the  peak 
computed  radial  stress  is  largest  with  the  cap  model,  and 
smallest  with  the  elas t i c-ideally  plastic  model.  This  result 
agrees  with  the  uniaxial  strain  stress  paths  shown  in  Fig.  (8). 
When  the  vertical  stresses  attenuate  to  less  than  200  psi,  the 
radial  stress  computed  with  the  elas tic-ideally  plastic  model 
is  the  largest,  a  result  again  suggested  by  Fig.  (8). 

Figure  (13)  shows  the  attenuation  of  the  peak  vertical 
and  horizontal  velocities  with  depth  at  the  same  500  psi 
pressure  contour.  As  was  the  case  with  the  vertical  stresses 
in  Fig.  (12),  the  peak  vertical  velocities  computed  with  the 
variable  moduli  and  elas tic-ideally  plastic  models  agree 
extremely  well,  reflecting  the  similarity  in  their  uniaxial 
strain  stress-strain  relations.  At  sufficient  depth  the 
additional  hysteresis  in  the  cap  model  produces  smaller 
vertical  velocities,  again  following  the  pattern  in  Fig.  (12). 
The  peak  horizontal  velocities  computed  with  the  three  models 
agree  fairly  well,  at  least  to  depths  of  z  -  200  ft.  It 
should  be  noted  that  the  horizontal  velocities  do  not  attenuate. 

The  attenuation  in  the  peak  vertical  and  horizontal 
displacements  at  the  500  psi  pressure  contour  is  shown  in 
Fig.  (14).  In  each  case,  the  results  are  fairly  close  and  the 
trends  are  similar  for  all  three  computations.  The  elastic- 
ideally  plastic  model  consistently  predicts  the  smallest 
horizontal  displacement.  It  should  again  be  noted  that  the 
horizontal  displacements  do  not  attenuate. 
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It  is  also  of  interest  to  compare  time  histories  of 
displacements,  velocities  and  stresses  at  points  in  the  medium. 
Figure  (15)  shows  vertical  displacements  versus  time  at  a 
near  surface  point  below  the  500  psi  surface  pressure  load 
(z  =  30  ft,  R  =  20/u  ft).  The  results  for  all  three  cases  are 
quite  comparable.  The  remarkable  agreement  between  the 
variable  moduli  and  elas ti c-ideally  plastic  models  may  again 
be  attributed  to  the  similarity  in  their  uniaxial  strain 
stress-strain  relations.  Figure  (16)  shows  horizontal 
displacement  versus  time  curves  at  the  same  point.  It  is 
seen  that  the  time  histories  for  the  variable  moduli  and  cap 
models  are  very  close  while  that  of  the  elas tic-ideally 
plastic  model  with  =  0.25  is  quite  similar  but  somewhat 

lower.  The  horizontal  displacement  versus  time  curve  for 
the  earlier  elas tic-ideally  plastic  model,  =  0.30,  is 

also  shown.  It  is  of  interest  to  note  the  considerable 
change  in  the  values  of  6^  when  is  changed  from  0.30  to 

0.25  in  the  e las  tic-idea lly  plastic  model.  What  is  significant 
here  is  not  the  change  in  the  value  of  ,  but  rather  the 

change  that  this  engenders  in  the  unloading  stress  path  in 
the  two  models.  In  Fig.  (8),  line  (l)  is  approximately  naifway 
between  line  (?)  (corresponding  in  =  0.30)  and  the  cap 

unloading  path.  Similarly,  the  horizontal  displacement 
computed  with  =  0.25  is  about  halfway  between  the  cap 

curve  and  the  =  0.30  curve.  Again,  the  importance  of 

modeling  the  correct  Stress  paths  in  both  loading  and  unloading 
in  uniaxial  strain  is  indicated. 


I 
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Figures  (17)  and  (18)  show  plots  of  vertical  velocity 
and  vertical  stress  versus  time  at  the  point  in  question.  In 
each  case,  the  comparable  curves  from  the  three  models  are 
very  similar,  in  fact,  they  are  almost  identical  at  early  times. 

The  time  history  of  radial  stress  computed  with  all 
three  models  is  shown  in  Fig.  (19).  Despite  the  numerical 
oscillations,  the  three  curves  are  seen  to  follow  distinct 
trends  after  the  first  few  oscillations.  The  cap  and  variable 
moduli  curves  are  close  to  one  another  and  significantly  above 
the  elas t ic-ideally  plastic  results  for  about  the  first  third 
of  the  elapsed  time  shown.  Thereafter,  the  radial  stress 
computed  with  the  variable  moduli  model  approaches  that  computed 
with  the  e las tic-ideally  plastic  model.  This  behavior  agrees 
with  that  shown  in  Fig.  (8)  where  the  variable  moduli  stress 
path  on  unloading  starts  close  to  the  cap  stress  path  and 
approaches,  but  is  always  to  the  right  of,  line  (T)  .  This 
constant  referral  to  the  stress  path  in  uniaxial  strain,  while 
extremely  useful,  is  not  the  whole  story  in  the  present  two 
dimensional  problem.  From  Fig.  (18)  it  is  seen  that  the 
vertical  stress  goes  essentially  to  zero  at  about  one  second 
while  the  radial  stress,  Fig.  (19),  continues  to  vary  at  much 
3  a ter  times . 

Finally,  the  he ri zontal  velocity  time  histories  at  the 
same  point,  computed  with  the  three  models,  are  shown  in 
Fig.  (20).  The  general  character  of  each  of  the  curves  is 
similar  and  the  cap  and  the  variable  moduli  curves  agree  most 
closely.  The  numerical  oscillations,  particularly  at  late  times. 


seem  to  be  smallest  in  the  variable  moduli  curve.  This  result 
is  anticipated  by  the  fact  that  the  variable  moduli  is  the 
only  model  of  the  three  which  contains  additional  hysteresis 
for  cyclic  loading  within  the  current  yield  surface. 

As  a  first  step  towards  extending,  the  present  study  to  a 
layered  geometry,  the  transmitting  bottom  boundary  in  the 
present  geometry,  Fig.  (1),  was  replaced  by  a  rigid  bottom 
boundary.  This  approximates  the  effect  of  an  underlying  layer 
of  hard  rock.  All  other  parameters  were  unchanged,  the 
elas tic-ideally  plastic  calculation  being  made  with  =  0.25, 

Preliminary  results  from  this  series  of  runs  are  shown  in 
Figs .  (21)  and  (22) . 

The  variation  of  peak  vertical  displacements  with  range 
and  overpressure  level  is  shown  for  the  rigid  bottom  calcu¬ 
lations  with  all  three  models  in  Fig.  (21).  The  previous 
results  from  Fig.  (10)  are  repeated  here  for  comparison. 

As  can  be  expected,  the  vertical  displacements  are  reduced 
by  the  presence  of  the  rigid  bottom  boundary.  In  addition, 
the  calculations  agree  extremely  well  with  each  other,  as  was 
true  of  the  earlier  transmitting  bottom  calculations. 

The  variation  of  peak  horizontal  displacements,  from  the 
rigid  bottom  calculations,  is  jhown  in  Fig.  (22),  super¬ 
imposed  upon  the  previous  results  of  Fig.  (11).  The  effect 
of  the  rigid  bottom  on  horizontal  displacements  is  not  as  easy 
to  predict  as  in  the  vertical  case.  The  horizontal  displacements 
are  increased  at  high  pressure  levels  (as  they  would  be  in  a 
fluid)  and  decreased  at  lower  pressure  levels  where  the  reflected 
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shear  wave  becomes  of  some  importance.  The  agreement  between 
the  results  predicted  with  the  various  models  is  nowhere  as 

good  as  it  was  previously  and,  in  fact,  there  again  appears 

to  be  a  difference  of  a  factor  of  two  between  the  elastic- 

ideally  plastic  model  calculation  (uj  th  =  0.25)  and  the 

variable  moduli  model.  The  cap  results  generally  fall  in 
between  the  other  two,  but  agree  much  more  closely  with  those 
from  the  variable  moduli  calculation.  The  analogy  to  the 
simple  uniaxial  strain  case,  which  was  so  useful  previously, 
is  less  satisfactory  in  explaining  the  current  results. 
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IV  CONCLUSIONS  AND  RECOMMENDATIONS . 

Some  interesting  conclusions  can  be  drawn  on  the  basis 
of  this  series  of  calculations  on  a  relatively  simple  geometry. 
First,  the  material  models  when  fitted  to  laboratory  data  gave 
comparable  vertical  results  (displacements,  velocities  and 
stresses)  in  each  case.  F  c  the  radial  quantities,  however, 
comparable  results  can  only  be  obtained  when  the  stress  paths 
in  a  loading-unloading  cycle  in  uniaxial  strain  are  also 
comparable . 

The  present  simulated  half-space  problem,  which  is  highly 

superseismic  (Mach  -  8  at  the  500  psi  level),  leads  to  states 

of  stress  which  are  very  close  to  those  in  uniaxial  strain. 

The  maximum  slope  of  the  peak  vertical  displacement  versus 

depth  curve,  Fig.  (14),  suggests  that  the  maximum  vertical 

strain  near  the  surface  (z  =  30  ft)  at  the  500  psi  overpressure 
* ) 

level  is  about  3%  ,  while  the  slope  at  the  corresponding  point 

for  the  peak  horizontal  displacement  versus  range  curve, 

*  *  ) 

Fig*  (11),  suggests  a  maximum  radial  strain  of  less  than  0.1% 

It  is  thus  seen  that  the  state  of  strain  is  nearly  uniaxial 
and  that  consequently,  any  model,  even  a  one  dimensional  model, 
which  closely  matched  the  experimental  uniaxial  strain  test 
stress-strain  curve,  Fig.  (2),  would  give  good  results  for 

It  should  be  noted  that  although  this  slope  is  a  strain¬ 
like  quantity,  it  is  not  a  true  strain  since  the  peak 
displacements  occur  at  different  times. 

The  maximum  hoop  strain  at  this  point  is  known,  i.e., 

6h 

£.=-#=  0.06%. 

<!>  R 


*) 

**) 
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vertical  stresses,  velocities  and  displacements.  A  one 
dimensional  model,,  however,  would  of  course  give  no  information 
on  the  small  (as  compared  to  the  vertical  effects)  but  real 
horizontal  velocities  and  displacements  which  o;ccur  in  this 
two  dimensional  prpblem.  » 

The  radial  stress,  velocity  and  displacement  quantities 
are  quite  sensitive  to  both  the  uniaxial  strain  test  stress- 
strain  curves  and,  in  particular,  to  the  stress  path  of  the 
material  in  a  load-unload  cycle.  The  radial  stress  time 

,  f  it 

history  is  largely  governed  by  the  vertical  stress  and  the 

i 

stress  path  of  the  material  in  uniaxial  strain.  The  radial 
stress  in  turn  affects  the  horizontal  velocity  and  hence^  the 
horizontal  displacement.  Consequently,  the  stress  path  of 
the  material  in  a  load-unload  cycle  plays  a  crucial  role  in 

1  i 

the  evaluation  of  the  horizontal  quantities, 

t 

i 

It  is  seen  from  the  results  that  the  horizontal  motion 
for. all  tb-ee  models  agreed  fairly  well,  although  the  results 

*  i 

for  the  elas t ic-ideally  plastic  model  were  smaller  than  those 

1  j 

for  the  others.  The' fact  that  the  three  models  agreed  as  well 
as  they  did  suggests  that  for  this  problem,  the  e las t ic-ideally 
plastic  model  would  be  a  reasonable  approximation  to  the’more 
sophisticated  models,  provided  that  the  stress  path  over  a 
uniaxial  strain  test  load- unload  cycle  were  also  modeled. 

•  f  , 

The  approximation  becomes  even  better  at  higher  stress  levels, 

i 

where  the  applied  stress  is  significantly  higher  than  thef 

>  ( 
deviatoric  stresses.  ’  f 
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The  results  point  out  the  importance  of  having  more  com¬ 
plete  laboratory  data  from  the  uniaxial  strain  laboratory  tests 
available  to  the  model  fitter.  While  the  vertical  effects  are 
primarily  governed  by  the  stress-strain  curve  for  hysteretic 
materials  of  this  type,  it  is  imperative  to  model  the  stress 
path  in  a  loading-unloading  cycle  for  a  good  definition  of  the 

i  i 

radial  effects.  This  means  that  in  addition  to  the  usual 
longitudinal  stress  versus  strain  data  from  a  uniaxial  test, 
lateral  stress  versus  longitudinal  stress  data  from  this  test 
is  also  urgently  required.  For  effects  in  the  superseismic 
range,  it  must  be  concluded  that  such  data  must  be  an  integral 
part  of  the  model  fitting  procedures. 

The  preliminary  results  for  the  case  of  a  rigid  bottom 
were  discussed  in  Section  III.  The  vertical  displacements 
agree  as  well  as  those  for  the  simpler  half-space  problem. 

■The  horizontal  displacements  on  the  other  hand  do  not  agree 
nearly  as  well  as  those  for  the  half-space  problem.  Once 
again,  the  elas t ic-ideally  plastic  model  gives  smaller 
horizontal  displacements,  sometimes  by  as  much  as  a  factor 
of  two.  Figure  (22)  serves  as(  a  warning  as  to  the  extent  of 
the  uncertainties  which  will  exist  in  the  prediction  of 
horizontal  motions  for  more  complicated  geometries. 

Ou  the  basis  of  the  present  study,  it  is  not  possible 
to  recommend  a  single  model  as  being  most  appropriate  for 
use  in  all  ground  shock  calculations.  For  the  simple  half¬ 
space  problem  considered,  all  three  models  gave  comparable 
results  as  long  as  the  unloading  stress  paths  were  somewhat 
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similar.  The  cap  model  is  the  only  one  which  satisfies  all 
theoretical  requirements  and  at  the  same  time  fits  all 
available  laboratory  data.  On  the  other  hand,  it  is  more 
difficult  to  fit  than  the  variable  moduli  model  which  also 
fits  all  the  laboratory  data.  The  variable  moduli  model, 
however,  is  theoretically  correct  only  for  cases  of 
proportional  loading  in  shear.  The  elas tic-ideal ly  plastic 
model  cannot  fit  all  of  the  laboratory  data.  The  use  of  a 
constant  Poisson's  ratio  and  a  nonassociated  flow  rule  in 
this  model  can  each  lead  to  theoretical  difficulties.  The 
computer  time  required  in  the  calculations  with  the  two 
advanced  models  (for  the  same  number  of  tJme  steps)  is 
similar  and  less  than  twice  that  required  for  the  elastic- 
ideally  plastic  model  computations. 

These  studies  are  currently  being  extended  to  layered 
geometries  and  other  materials.  In  addition,  effects  in  the 
transseismic  and  possibly  subseismic  regions  will  also  be 
considered,  in  an  effort  to  arrive  at  specific  model 
recommendations  for  a  wide  variety  of  problems  of  interest. 


TABLE  I  -  Summary  -  Advantages  and  Disadvantages  of  Each  Model. 
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FIG.  8  COMPARISON  OF  CAP  AND  ELASTIC  -  IDEALLY  PLASTIC  MODELS 
STRESS  PATH  FOR  UNIAXIAL  STRAIN  TEST 
LOADING  TO  400  PS  I  AND  UNLOADING. 
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FIG.  9  COMPARISON  OF  CAP  AND  ELASTIC  -  IDEALLY  PLASTIC  MODELS 
STRESS  PATH  FOR  UNIAXIAL  STRAIN  TEST 
LOADING  TO  IOOO  PS  I  AND  UNLOADING. 
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